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ABSTRACT 

Starting from the observed distribution of galaxy clusters in redshift space we 
use a two-step procedure to recover their distances and peculiar velocities. Af- 
ter statistically correcting for the unobserved cluster distribution in the zone 
of avoidance < 20°) and also for a smooth absorption at higher |6|'s, we 
use a dynamical iterative algorithm to recover the real-space cluster positions 
by minimizing the redshift space distortions. The whole procedure assumes 
that clusters trace the mass, that peculiar velocities are caused by gravity and 
that linear perturbation theory applies. The amplitude of the cluster dipole 
measured in the 3D space turns out to be ~ 23% less than that measured in 
redshift space. In both cases the dipole direction is aligned with the Cosmic 
Microwave Background dipole within ~ 10°, taking into account the Virgocen- 
tric infall component of the Local Group [LG hereafter] motion. Observational 
errors, limitations in the reconstruction procedure and the intrinsic cosmolog- 
ical variance, which is the dominant source of uncertainty, render a stringent 
determination of the /3(= f2°'^/6) parameter whose central value turns out to be 
(3 0.2 while its total uncertainty is ±0.1. This implies that for a cluster-mass 
bias parameter of ~ 5, a flat Universe is not excluded, contrary to previous 
cluster-dipole z-space analysis. A more stringent determination of (3 will be 
obtained from the analysis of the peculiar velocity field in a forthcoming paper. 

radiation sea, towards I = 277° and h = 30° (Lubin &; 
Villela 1986; Kogut et al. 1993) 

If linear gravity is responsible for the observed pe- 
culiar velocities then according to the linear instability 
theory (of. Peebles 1980) the peculiar gravitational ac- 
celeration should be aligned and proportional to the pe- 
culiar velocity; the constant of proportionality being a 
measure of the present-time growth rate of mass fluc- 
tuations and therefore a measure of the cosmological 
density parameter, The gravitational acceleration 
of any galaxy (or other extragalactic object) can be es- 
timated by calculating the dipole moment of the distri- 
bution of mass surrounding it. Since, however, the latter 
is unknown one has to resort in estimating the dipole 
moment of the distribution of luminous extragalactic 
objects using some simplifying assumption about the 
relation between fluctuations in the matter and light 
distributions; which is non other than the linear-biasing 
assumption (cf. Kaiser 1984). In this picture the galaxy 



1 INTRODUCTION 

The study of cosmic peculiar velocities is a relatively re- 
cent and very exciting field that could give great insight 
to the origin of the large-scale structures in the Uni- 
verse. With the advent of new observational techniques 
and telescope power it has become possible to relate 
distance dependent with distance independent quanti- 
ties and thus determine the relative distance for a large 
number of galaxies; the most frequently of these being 
used are the Tully-Fisher relation for spirals and the 
Dn — a relation for elliptical galaxies (cf. Burstein 1990 
and Dekel 1994 for comprehensive reviews). 

The best evidence that cosmic structures can have 
significant velocities, above their cosmological expan- 
sion velocities, came from the interpretation of the CMB 
temperature dipole anisotropy as a Doppler effect, orig- 
inating from the motion, with vlg = 622 ± 20 km/sec, 
of the Local Group of galaxies in the isotropic CMB 
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fluctuations, 6g, are related to those of the matter dis- 
tribution (5 by a constant factor, the biasing parameter, 

6g = b6. 

Further difficulties arise from the fact that one can 
observe galaxies or other extragalactic objects limited 
in magnitude or flux, which then implies that distant 
contributions to the dipole could be missed if the char- 
acteristic depth of the galaxy catalogue is less than the 
convergence depth of the dipole. Moreover, most cat- 
alogues of extragalactic objects have also limited sky 
coverage which is usually due to light absorption near 
the Galactic plane. 

Up to now the dipole of various populations of ex- 
tragalactic objects has been determined: optical galax- 
ies (Lahav 1987; Plionis 1988; Lahav, Rowan-Robinson 
&; Lynden-Bell 1988; Lynden-Bell, Lahav &; Burstein 
1989; Hudson 1993b), IRAS galaxies (Meiksin fc Davis 
1986; Yahil, Walker fc Rowan-Robinson 1986; Villum- 
sen &; Strauss 1987; Strauss &; Davis 1988; Strauss et 
al. 1992; Rowan-Robinson et al. 1990; Plionis, Coles 
&; Catelan 1993), X-ray active galactic nuclei (Miyaji 
fc Boldt 1990); X-ray clusters (Lahav et al. 1989) and 
Abell clusters (Scaramella, Vettolani &; Zamorani 1991 
[hereafter SVZ91]; Plionis &; Valdarnini 1991 [hereafter 
PV91]). In all cases the dipole moment was found to 
be quite well aligned with the CMB dipole suggesting 
that gravity is indeed responsible for the Local Group 
motion and that light traces mass. 

An interesting historical fact, originating from these 
studies, is that the estimated dipole convergence depth 
has been a function of cosmic volume sampled; the 
deeper the catalogue the larger the dipole convergence 
depth. This is true out to the largest depths, traced by 
the Abell/ AGO cluster catalogue, for which Rconv ~ 
17000 km/sec (PV91 and SVZ91). This implies that the 
apparent dipole convergence of the shallower galaxy cat- 
alogues (being optical or IRAS) is probably artificial, 
determined by their limiting depth. In fact a careful 
investigation of the QDOT-IRAS sample showed some 
evidence for a contribution to the dipole from scales 
comparable to the cluster Rconv (Plionis, Goles &; Gate- 
Ian 1993). 

A further interesting outcome of the Abell/ AGO 
cluster dipole analysis is the fact that there seems to be 
a coherent d&poVe anisotropy out to ~ 17000 km/sec; i.e., 
the differential gravitational acceleration induced on the 
LG by the distribution of clusters in large equal- volume 
shells is roughly aligned with the GMB dipole in each 
shell (PV91). Furthermore, there is strong evidence for 
the existence of coherent large-scale galaxy fiows in the 
local Universe, extending from the Perseus-Pisces region 
on the one side to the Hydra-Gentaurus/Great Attrac- 
tor region on the other. These fiows are well established 
both for elliptical (cf. Lynden-Bell et al. 1988) and spiral 
galaxies (cf. Rubin et al. 1976; Dressier &; Faber 1990; 
Mould et al. 1991; Willick 1990; Han fc Mould 1990; 



Gourteau et al. 1993; Hudson 1994; Mathewson &; Ford 
1994) and point in the general direction of the GMB 
dipole. These results put together present a consistent 
picture in which the LG participates in a large-scale 
bulk motion induced by gravity, encompassing a large 
volume of radius ~ 6000 — 15000 km/sec (see the review 
of Dekel 1994). Although the extent of this volume may 
still be under debate the existence of the bulk fiow seems 
to be certain. 

This picture has been recently challenged by Lauer 
&; Postman (1994) [hereafter LP94] who have extended 
the cosmic fiow studies to very large scales, those traced 
by galaxy clusters, and who find that the LG motion 
with respect to the frame defined by Abell/ AGO clus- 
ters, within 15000 km/sec, moves in a direction very 
different than that of the GMB dipole which then im- 
plies that, if the GMB dipole is a Doppler effect, the 
whole cluster frame is moving with respect to the GMB 
rest-frame with ~ 700 km/sec. Such a large velocity on 
such large scales has put into despair structure forma- 
tion model-builders, since no model can easily accom- 
modate such velocities on large scales (cf. Strauss et al. 
1994; Feldman fc Watkins 1994; Jaffe fc Kaiser 1994). 

Furthermore and more importantly it puts into 
doubt the gravitational instability picture because on 
the scales traced by clusters of galaxies linear the- 
ory should apply, which predicts that acceleration and 
velocity are aligned and as discussed previously, the 
gravitational acceleration of the LG of galaxies, de- 
termined from the cluster distribution in z-space and 
within ~ 25000 km/sec, is found to be very well aligned 
{SOcmb < 15°) with the GMB dipole (PV91; SVZ91). 

So we are left with the following puzzling pic- 
ture: the acceleration determined from the distribution 
of galaxies (within ~ 10000 km/sec) and of clusters 
(within ~ 25000 km/sec) is well aligned with the GMB 
dipole (which is in itself a strong indication that the 
GMB dipole is Doppler generated) while the LG veloc- 
ity (as determined by LP94), with respect to the cluster 
frame within 15000 km/sec, points to a direction al- 
most perpendicular to the GMB dipole direction. One 
would then be forced to explain the alignment of the 
LG gravitational acceleration, determined from the dis- 
tribution of clusters and galaxies, with the GMB dipole 
direction as a 'cosmic' coincidence. However, the ran- 
dom joint probability of having the cluster dipole, the 
optical and IRAS galaxy dipoles aligned with the GMB 
dipole within a few degrees (~ 10°, ~ 10° and ~ 25° 
respectively), assuming that they are independent, is 
< 3 X 10~^. Furthermore, as we discussed previously, 
PV91 found evidence for a coherent dipole anisotropy, 
in which the differential cluster dipole in each equal- 
volume shell roughly points towards the GMB dipole. In 
fact, the dipole in the first shell (ii^ 100 h~^ Mpc) has 
89< 20° while in the most distant shell (140< R< 160 
h~^ Mpc), for which the dipole is aligned with the GMB, 
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it has 69^ 10° and the distance between the centers of 
the two shells is ~ 90 Mpc, a distance at which 

^cc('') <C 1, where ^cc is the spatial cluster-cluster cor- 
relation function (note that even the edges of the two 
shells are ~ 43 Mpc away). Therefore, these two 

volumes could be considered independent which would 
then decrease even further the probability of random 
alignment with the CMB dipole to < 8 x 10"*! 

We attempt to throw light to this paradox by using 
linear perturbation theory and a dynamical reconstruc- 
tion algorithm (similar to that of Strauss &; Davis 1988) 
to determine the 3D cosmic density field, traced by the 
Abell/ACO galaxy clusters within < 20000 km/sec, 
with the aim of investigating 

(i) whether the large cluster-dipole convergence depth, 
its asymptotic amplitude and the alignment of the 
z-space cluster dipole with that of the CMB, are 
artifacts of redshift space distortions, 
(ii) the peculiar velocity field, within 20000 km/sec, 
predicted by linear theory and comparing it directly 
with that derived by LP94. 

In this first paper we will address the first set of 
questions, outlined above; while the second will be ad- 
dressed in a forthcoming paper. The plan of this paper 
is as follows: In section 2 we present an extended dis- 
cussion of the two subsamples used, of their biases and 
homogenization procedure. In section 3 we present both 
the zone of avoidance filling procedure, the dynamical 
real-space reconstruction algorithm and its reliability 
tests. In section 4 we present the cluster dipole esti- 
mate, possible systematic effects and an extended error 
analysis. In section 5 we derive interesting cosmologi- 
cal parameters from the dipole analysis and finally in 
section 6 we present our conclusions. 

2 THE CLUSTER SAMPLE 

Our present analysis is based on a cluster sample ex- 
tracted from the original Abell sample (Abell 1958) and 
from its southern extension, the AGO cluster catalogue 
(Abell, Corwin &; Olowin 1989). The north declination 
limit of the AGO sample is 6 = —17, so that the two 
catalogues overlap in the strip —27 < (5 < —17. The 
overlapping region has been used by Abell, Corwin &; 
Olowin (1989) to calibrate the J magnitudes of the 
AGO to the Abell magnitude system. However, system- 
atic differences in the two calibrated catalogues have 
been noticed by AGO (1989), Batusky et al. (1989), 
Scaramella et al. (1990) and PV91. The problem of ob- 
taining a statistically homogeneous catalogue from the 
composite Abell/ AGO sample, is particularly important 
in dipole studies since artificial systematic differences 
(density gradients or global density variations) can en- 
hance or even produce a dipole signal. Although this 
issue has been addressed by a number of authors (cf. 



Scaramella et al. 1990 and PV91) we will also discuss 
it, following the PV91 line of reasoning. 

Gluster distances, r, are estimated from redshifts 
using the standard relation: 

^ = H./il+z) + - {' - V^^^)] '(1) 

where c is the speed of the light and qo the deceleration 
parameter. The number of clusters within our subsam- 
ple is not fixed since the distances slightly depend on 
Qo which we do not fix a priori. We also merge cluster 
pairs with relative distance < 3 Mpc into unique 
'clusters' having as position their center of mass and as 
mass their combined total mass. We made this approx- 
imations to avoid nonlinear effects on small scales for 
which our reconstruction method fails to account for. 
The combined subsample we consider is composed by: 

(i) Abell clusters in the roughly volume-limited region 
r < 250 h-^ Mpc, |6| > 13° and with mio < 17. 
mio is the magnitude of the tenth brightest galaxy 
in the cluster in the magnitude system corrected 
according to PV91. This sample is 100 % red- 
shift complete. Note that the number of objects 
in this subsample varies between 269 and 276 for 
Qo = 0.2 and 0.5, respectively. Out of these we have 
merged 5 cluster pairs, 
(ii) AGO clusters with r < 250 h'^ Mpc, |6| > 13° 
and mio < 17. Among the 217 objects within the 
sample ~ 77% have measured redshift. For the re- 
maining objects the redshift is determined using the 
mio — z relation used by PV91. Note that we have 
merged 7 nearby cluster pairs. 

In order to obtain a statistically homogeneous all- 
sky sample of clusters we quantify separately the obser- 
vational biases, present in the Abell and AGO subsam- 
ples, and then we calibrate these subsamples applying 
a homogenization procedure. 

2.1 The Selection Biases 

At low galactic latitudes, absorption of optical photons 
makes it unlikely to observe luminous objects. In the 
galactic strip |6| < 20°, to which we refer as the Zone 
of Avoidance (ZoA), there are only 20 clusters in our 
subsample. The galactic absorption outside the ZoA has 
been always found to be consistent with a cosecant law: 

log P{b) = a{l - cosec|6|), (2) 

where P{b) gives the probability that a cluster in the 
range 6 to 6 + A6 would have been included in the cat- 
alogue. There is no general agreement on the precise 
value of a. Bahcall &; Soneira (1983) and Postman et 
al. (1989) found 0.3 for Abell clusters, Batuski et al. 
(1989) found a ~ 0.2 for AGO clusters while LP94 
found an even lower value (a ~ 0.15) for both AGO 
and Abell clusters contained in a volume-limited 15000 
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km/sec sphere. To allow for this experimental uncer- 
tainty we will consider, in what follows, two different 
sets of values: (a,a)= (0.3,0.2), that better matches the 
behaviour of the number density of our cluster as a func- 
tion of the galactic latitude and (0.2,0.2). The first value 
in each set refers to the Abell sample while the second 
to the AGO one. 

The finite depth of the cluster sample was modeled 
with the redshift selection function that was estimated 
as in Postman et al. (1989). We modeled the probability 
that a cluster in the range z to z + Az would have been 
included in the cluster catalogue as: 



that the two catalogues are simply radially inhomoge- 
neous, we computed the following weighting function: 



P{z) 



{1 



(3) 



i{ z < Zc 
exp[—z/zo) i{ z > Zc 

where Zc is the redshift up to which the space density of 
clusters remains constant (volume-limited regime). We 
obtain Zc ~ 0.0787 and 0.0664 for the Abell and AGO 
subsamples respectively. 

No appreciable declination dependence has been 
detected in the Abell/ AGO subsample (but see 
Scaramella et al. 1990 for a possible zenithal depen- 
dence). 

2.2 Homogenization of the Abell/ AGO 
sample 

In this work we will follow the PV91 homogenization 
scheme, used to minimize the density variations of the 
Abell and AGO subsamples. The estimated observa- 
tional biases allow us to assign the following weight to 
clusters: 

1 



Ui{z,b) 



(4) 



P{zi)P{bi)- 

The individual number density of clusters of the two 
subsamples is therefore n = ^^.Wi[z,b), where the 
sum extends over all clusters within the sampled volume 
V. Since at low galactic latitudes the patchy galactic 
absorption is not well described by eq.(2) we compute 
n above |6| = 30°. We divided V into 15 equal vol- 
ume intervals in which we estimated the Abell and AGO 
cluster number density separately. Due to the complex 
geometric boundaries, the volumes were computed by 
Montecarlo calculations. 

Within ~ 200 Mpc the space density of AGO 
and Abell is roughly constant with n ~ 2.5 ± 0.3 x 10~^ 
and 1.5±0.2xlO"'^ /i^ Mpc"^ , respectively. Errors were 
computed from the density fluctuations in the different 
bins. Although beyond ~ 200 Mpc the density er- 
rors of the Abell sample remain roughly the same, those 
of the AGO increase dramatically (~ lO"'^ /i^ Mpc"^) 
implying large density fluctuations due to undersam- 
pling. 

An approximate homogenization of the two samples 
was then done by equating the number density of the 
two catalogues in each equal volume region. Assuming 
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Abell) 



i{6> 
if <5 < 



-17 
-17 



(5) 



,ACO) 

where Ti is the distance of the i"" bin. No significant 
difference was found when using the two different sets 
of absorption coefficients. Note, however, that this is 
not the only possible homogenization scheme since we 
do not know a priori which is the true value of the 
cluster density. We could, therefore, homogenize the two 
samples by using the inverse of eq.(5), yVrei{r)~^ . As we 
will show in section 4 our main results remain stable on 
the choice of the homogenization scheme. 



3 THE RECONSTRUCTION METHOD 

The whole reconstruction procedure can be schemati- 
cally represented by the following diagram: 



Zofes 



Zid 



3D„ 



Starting from the observed sky redshift space distribu- 
tion of galaxy cluster, Zobs, we generate a synthetic clus- 
ter distribution to account for unobserved clusters. This 
procedure allows us to reconstruct the whole-sky clus- 
ter distribution in the redshift space, Zid- An iterative 
reconstruction algorithm, similar to that proposed by 
Strauss &; Davis (1988), is then applied to minimize the 
redshift space distortion, allowing us to recover the real- 
space distribution of the observed clusters 3Do6s and 
therefore also their peculiar velocities. The intrinsic re- 
liability of the entire procedure has been tested using a 
simulated catalogue of clusters kindly provided to us by 
S.Borgani. 

3.1 z-Space reconstruction scheme 

To reconstruct the whole-sky redshift space cluster 
distribution, Zid, we adopted a phenomenological ap- 
proach. The basic idea is to fill the artificially (due to 
Galactic absorption) low density regions of the surveyed 
volume with a synthetic cluster population having the 
same clustering properties as the distribution of the real 
clusters and then generate many Monte-Garlo realiza- 
tions of this population. To this end we divide our pro- 
cedure into two steps; one to reconstruct the cluster 
density field above |6| = 20° and one to fill the ZoA. 

3.1.1 \b\ > 20° 

The first step is to recover the cluster distribution for 
|6| > 20°. We divided our volume into two regions: an 
inner sphere with a radius of Rj„=200 Mpc, which 
we consider as the region of reliable determination of the 
cluster density field, and an external region Rj„ < ?• < 
250 h-^ Mpc. 
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In the inner part, in which P{z) ~ 1, we used a 
Montecarlo rejection method to generate a population of 
synthetic clusters distributed according to the adopted 
P{b) probability function and weighted them by 1/P[z) 
with the further constraint of being spatially clustered 
around the real clusters according to the observed spa- 
tial cluster-cluster correlation function (cf. Bahcall &; 
Soneira 1983; Postman et al. 1992; Plionis, Valdarnini 
&; Jing 1992). The following weights are then assigned 
both to real and mock clusters: 



Finally, the total number of synthetic clusters outside 
the ZoA is determined by requiring the overall number 
density of clusters to be: 

n = ^'^Wi{z,b), (7) 

i 

where Wi[z,b) = Wrei{ri)wi[z,b). This requirement is 
fulfilled by imposing: 

real real synth 

where the first two sums extend over real clusters and 
the last one over synthetic objects. Beyond 200 
Mpc the radial selection function, mainly for the AGO 
sample, falls below one roughly exponentially. Therefore 
weighting clusters that sparsely trace the density field 
with 1/P[z) can introduce non negligible shot noise er- 
rors. To avoid this problem and keeping in mind that we 
will limit our analysis within Rj„ we reconstructed the 
cluster distribution in the external regions as following: 
the mock cluster distribution was Montecarlo generated 
to account for the P{b) selection with no spatial corre- 
lation with real clusters; then the cluster redshift was 
determined according to P{z) but, instead of weighting 
them by 1/P[z), we generated 1/P[z) synthetic clusters 
at the same b and z but with a randomly chosen I. Both 
real and synthetic clusters in this external region were 
simply weighted by Wrei{r) and their total number was 
determined by equating the local density with that of 
the central region. 

3.1.2 Filling the ZoA 

The second step is to fill the ZoA. The method we used 
is close to that proposed by S. Faber and used by Yahil 
et al. (1991). We divided the equatorial strip —20° < 
b < +20° into 18 bins of 20° in longitude. In each bin 
we divided the distance range into bins of 2000 km/sec. 
The clusters were then sampled in the two northern and 
southern adjacent strips, defined so that their joint so^A 
angle is equal to that of the ZoA. This requirement leads 
to the following limiting latitude for the strips: 

\hii^\ = 90° - cos"^ [2 X cos 70°] = 43.16° . 



The number of synthetic clusters in each ZoA bin was 
then set equal to the sum of the clusters (real and syn- 
thetic) found in the two bins of the northern and south- 
ern, adjacent to the ZoA, strips; while their position, 
inside the ZoA bin, was randomly assigned. Finally, the 
real clusters in the range 13° < |6| < 20° are inserted 
and the total number of synthetic ZoA clusters is ad- 
justed by subtracting a few at random until the average 
number density of real clusters is reached. 

3.1.3 The Zid Cluster Distribution 

In what follows we will refer to the above 2-step 
^-reconstruction scheme as the Randomized Standard 
Cloned Mask; [RSCM]. At the end of the whole pro- 
cedure nearly 50 % of the cluster population is com- 
posed by synthetic objects. A typical Zid reconstruction 
is shown in panels (a) and (c) of Figure 1. Black dots 
represent real clusters while open dots are synthetic 
objects. The figure displays the projection of the re- 
constructed cluster distributions onto the supergalactic 
plane (Xs„p-Ys„p projection; panel a) and orthogonal to 
it (Xs„p-Zs„p projection; panel c). A simple visual in- 
spection reveals that the synthetic clusters fill the ZoA 
and smoothly reproduce the clustering observed in the 
two nearby galactic-latitude strips. Outside the ZoA, as 
expected, synthetic objects are clustered with the true 
ones and their density decreases towards the galactic 
poles. 

However, since each cluster has a different weight a 
more relevant visualization of the z-space cluster distri- 
bution is presented as smoothed density maps in panels 
(6) and (d) where we have overlayed a 20 x 20 grid (each 
cell being 2000 km/sec wide) on the cluster distribution 
of panels (a) and (c) and smoothed the projected distri- 
bution using a 2D-Gaussian kernel with Rsm = 1 cell, 
while weighting each cluster by Wi. What, in fact, we 
present in these figures is the 'mean' smoothed density 
field; an average over 10 random realizations of the z- 
space reconstruction procedure. Note that in panels (6) 
and (d) the cell units used have and 20 corresponding 
to the Cartesian supergalactic coordinate -20000 and 
20000 (in units of km/sec) respectively. It is interesting 
to note in panel (d) [Xs„p-Zs„p projection], the exis- 
tence of a 'cross'-like structure with low-density cylin- 
drical regions, passing through the cluster distribution 
and having a length of ~ 20000 km/sec. A more detailed 
discussion of the apparent structures will be presented 
in section 3.4. 

3.1.4 Variants of 'R.'&CM 

For the ZoA filling method we can use two variants of 
the previously described method. In the first variant, 
which we call Randomized Heavy Cloned Mask [RHCM], 
we fill independently the northern and southern parts of 
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the ZoA by applying the [RSCM] to the northern and 
southern adjacent galactic strips (20° < |6| < 43.16°) 
separately. With this mask the randomization applies 
to smaller bins and the cloning effect is expected to be 
enhanced. 

In the second variant, which we call [Randomized 
Thin Cloned Mask [RTCMJ), we identify the ZoA with 
the equatorial strip |6| < 13°. Then we consider the 
two adjacent strips defined so that their individual so^A 
angles are equal to that of the ZoA. The [RSCM] is then 
applied by half counting the number of clusters in the 
adjacent bins. The purpose of this scheme is to account 
for the spatial clustering on larger volumes. 

Finally, we also use an alternative scheme to re- 
construct the r > 200 Mpc region which is based 
in extending the reconstruction technique used for the 
r <200 volume to the external region. The ZoA 

filling method used is the standard one, described in 
the previous subsection. We call this scheme Modified 
Standard Scheme [MSS]. 

3.2 3D— Space Reconstruction Scheme 

Our 3D reconstruction scheme is based on the assump- 
tion that the peculiar velocities of clusters are caused 
by gravitational instability and that linear perturbation 
theory applies on the scales relevant to the cluster dis- 
tribution. The distribution of clusters in redshift space 
differs from the true three-dimensional one by a nonlin- 
ear term in the redshift-distance relationship: 



In order to apply linear theory we need to smooth 
the discrete cluster density field on an appropriate scale, 
where non- linear effects could be important. Since we 
are using galaxy clusters to probe the density and the 
peculiar velocity fields, a natural smoothing scale, large 
enough for linear theory to be valid, would be the clus- 
ter correlation length (i.e. the distance To at which the 
spatial cluster-cluster correlation function is unity). Re- 
cently Croft &; Efstathiou (1994), analysing large N- 
body simulations, found that the velocity fields traced 
by galaxy clusters are highly non-linear below ~ 10 h~^ 
Mpc. Therefore we allow the smoothing radius to vary in 
the range [10, 30] h~^ Mpc. Although the cluster pecu- 
liar velocity field, the analysis of which will be presented 
in a forthcoming paper, does depend on the choice of the 
smoothing length, the dipole parameters depend only 
weakly on its choice, as we will see below. 

The discrete cluster distribution was smoothed with 
a top hat window function which is equivalent to smooth 
the peculiar acceleration with: 



W-(|r-r|) 



^ if|r-r|<.3 
1 if Ir — rl > Ts 



(12) 



where Ts is the smoothing length. 

Since Zid is reconstructed up to Rmax = 2hOh~^ 
Mpc, the integral in eq.(lO) is converted into a sum over 
all the observed and synthetic clusters within Rmax'- 



v(r) 



(13) 



cz = ifo|r| + [v(r)-v(0)] • — , 



(9) 



^(^)^(,-,).(.«+4 



where r is the position of the generic cluster, v(r) its pe- 
culiar velocity, v(0) the peculiar velocity of the observer 
placed at the centre of the coordinate system and Ho is 
the Hubble parameter. 

3.2.1 Linear Theory 

In linear theory the peculiar velocities are proportional 
to the peculiar acceleration (cf. Peebles 1980) 

Hof{Qo) I r 



.(r) 



47r 



—rrd r 



(10) 



where 6{i) = [p{i) — Pb]/ Pb is the mass density fluctu- 
ation about the mean pb, /(fio) ~ n^ *^ and no is the 
cosmological density parameter at the present epoch. 
We assume that cluster density fluctuations are related 
to the mass fluctuations by a constant linear biasing 
factor: 



6,{i) = bj{i), 



(11) 



where be is the bias parameter of clusters with re- 
spect to the mass. We can therefore replace 6{i) with 
(5c(r) in eq.(lO) provided that we substitute f{Oo) with 

nr/6c (=/3) 



The average cluster density Uc is obtained by divid- 
ing the average cluster mass density pc with Mc (mass 
of Coma cluster); pc is given by: 



(14) 



where Mi is the mass of the generic cluster. As in PV91 
and SVZ91 we assume that cluster masses are propor- 
tional to the the Abell-catalogue listed number of galax- 
ies per cluster. Mi oc NA,i (and Nc = 106); while the 
mass of the synthetic clusters was set equal to the aver- 
age mass of the real clusters. The last term corrects for 
the net gravity of the homogeneous density background. 
The correction is necessary since the summation calcu- 
lates the contribution from the density rather than from 
overdensity. Eq.(13) measures the true peculiar velocity 
only if density inhomogeneities outside the surveyed vol- 
ume can be ignored. 

As shown by Vittorio &; Juskiewicz (1987) and 
Juskiewicz, Vittorio &; Wyse (1990) in a multipole ex- 
pansion of the external gravitational field the leading 
term is the dipole. Its effect is naturally removed if the 
peculiar velocities are evaluated in the LG frame be- 
cause the dipole contribution is constant for all objects 
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and does not affect the relative peculiar velocities. The 
only disadvantage of the LG frame is that the correction 
for the solar motion relative to its barycenter is uncer- 
tain by ~ 50 - 100 km/sec (Yahil et al. 1977) which, as 
we will see, is smaller than the errors on the predicted 
LG velocity. For this reason we used eq.(13) to evaluate 
the cluster peculiar velocity in the LG frame. The higher 
order terms in the multipole expansion of the external 
gravitational field fall off' faster than the dipole with dis- 
tance and are generally negligible for distances smaller 
than RmaxI'i (Yahil et al. 1991). Possible systematic 
errors when computing peculiar velocities are therefore 
expected only at larger distances and should not ap- 
preciably affect the reconstruction reliability within 200 
h-^ Mpc. 



3.2.2 The Reconstruction Algorithm 

To reconstruct the 3D density field we used an iterative 
algorithm (which we call Linear Iterative Reconstruction 
Algorithm, /i/iJM/ hereafter), which is close to those de- 
scribed by Strauss fc Davis (1988), Yahil et al. (1991) 
and Hudson (1993a). The idea is to start from the ob- 
served redshift space cluster distribution and attempt to 
recover the cluster true positions by iteratively comput- 
ing their peculiar velocities according to eq.(lO). This 
method is self consistent as long as the density field we 
are dealing with is linear, the clusters trace the mass and 
the fluctuations responsible for the peculiar motions of 
objects are within the volume sampled. Note also that in 
order to avoid problems related to the definition of dis- 
tance and to the uncertainties on the Hubble constant, 
cluster distances are in km/sec units. 

The steps of our iterative scheme are the following: 

(i) All the clusters are initially placed at their observed 
distance, rl"'' = cz, with no peculiar velocity and 
with an arbitrary value for the input /3 parameter. 
The index refers to the zeroth iteration. 

(ii) The weighting function ^(rj) is computed and its 
value is kept constant in the subsequent iterations 
(i.e. the selection function is not upgraded). This 
speeds up the algorithm while it also does not lead 
to appreciable errors, since P{^z) K, 1 within 20000 
km/sec. 

(iii) The average density Uc within the sampled volume 
is computed and its value is updated. 

(iv) The window function VF(|rj'''' — r'*'!) is computed 
for each object. 

(v) The peculiar velocity of the i-th cluster v'''''(rj) is 
computed according to eq.(13) taking into account 
both, the clusters within Rmax and those that dur- 
ing the iterations were placed beyond the volume 
limits (we will refer to these as 'passive' objects). 

(vi) Radial positions of 'active' clusters are updated. 
Their new positions are given by 



^<'=+^'=--[v(r)W-v(0)W].^. (15) 

Positions of 'passive' objects are not updated, 
(vii) For each real active cluster we compute the dif- 
ference of its position between the last two itera- 
tions. If the difference is larger than 0.5% r, then 
we jump back to step (iii). The convergence is typ- 
ically reached within 10 iterations. 

Since the smoothing scheme treats clusters as 
spheres of radius Ts, clusters in the spherical shell 
|r — Rmax\ < need to be weighted properly by com- 
puting how much of the sphere lies within Rmax- But 
since we will estimate the dipole and the cluster velocity 
field only within 20000 km/sec and since we regard the 
cluster distribution beyond this radius as an improve- 
ment over a simple homogeneous distribution, we do not 
need to implement this correction. 

Our iterative reconstruction algorithm can fail to 
recover the true 3D object position within high clus- 
ter density regions ('triple value regions') (cf. Yahil et 
al. 1991; Hudson 1993a). In our case, however, the ob- 
jects are so sparsely distributed (which is the reason 
why the smoothing scale used is relatively large) that 
this problem is negligible. Obviously, the sparseness of 
the clusters introduces different problems, as we will 
discuss below. 

3.3 Testing the reconstruction reliability 

To test the intrinsic reliability of the whole reconstruc- 
tion procedure and to optimize its performance we con- 
sidered a catalogue of mock clusters extracted from a 
simulation performed by Borgani et al. (1994) in which 
the present time cluster positions were determined by 
the Zel'dovich approximation in a standard CDM model 
Universe. The mock clusters were identified as peaks 
above a density threshold which was chosen so that 
their space density matches that of the observed clusters 
(~ 2 X 10~^h^ Mpc~^). We considered a spherical vol- 
ume centered around a LG-like observer (see Borgani et 
al. 1994). Each cluster peculiar velocity was determined 
according to linear theory after having computed the 
gravitational acceleration generated by all the clusters 
within the sampled volume. This ideal cluster distribu- 
tion was then degraded by Montecarlo rejecting clusters 
according to the probability distributions P{^z) and P{1>) 
and by devoiding the ZoA. Finally, an artificial Gaus- 
sian distributed noise was added to the cluster velocity 
to mimic experimental errors in their measured redshift. 
According to observational indications we used a vari- 
ance of 200 km/sec (Strauss, Cen &; Ostriker 1994). Fi- 
nally we applied the RTCM together with the LIRM to 
recover the original object distribution and velocities. 

Three different indicators were used to assess the 
intrinsic reliability of the reconstruction procedure: 
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• The difference between reconstructed and true pe- 
culiar acceleration acting on the observer. 

• The volume averaged discrepancy between the true 
and the reconstructed position: (ej) = {\itrue,i — 

^rec.,i\/\^true,i\^ 

• The spatial cluster-cluster correlation function. 

The reconstructed peculiar acceleration vector is 
very close to the original one. The per cent discrep- 
ancy between the amplitudes of the true and the re- 
constructed peculiar velocities, based on 10 different re- 
constructions, was 3 ± 6%. The two vectors were mis- 
aligned by 2° ± 5°. In Figure 2 we present the distri- 
bution of the ei values which peaks at ^ 1% and has a 
mean value of (ej) ~ 1.6% (it is skewed towards large 
errors); while the probability of having > 3% is quite 
low. We found, as expected, that depends on the ini- 
tial displacement of the i-th cluster with respect to its 
true position \itrue,i — i^°^\ - The reliability of the method 
greatly improves if clusters are initially placed close to 
their real positions, which would be possible if reliable 
redshift independent distance measurements were avail- 
able for some nearby clusters. Finally, the spatial two 
point correlation function of the reconstructed cluster 
distribution turns out to be indistinguishable from the 
true one. 

We stress that these tests are aimed at assessing the 
intrinsic reliability of the reconstruction method; clus- 
ter velocities are computed according to linear theory 
and no allowance for contribution from external gravi- 
tational fields is made (a set-up that needs not be true 
in the real Universe). To get an estimate of the actual 
uncertainties in the whole procedure we need to account 
for several other effects, which will be quantified in Sec- 
tion 5. The intrinsic uncertainty will be assumed to be 
independent from the other source of errors and will be 
added in quadrature to compute the total uncertainty 
on the peculiar velocities. 

3.4 Smoothed cluster density maps in z- 
and 3D— space 

In Figure 3a and b we plot the reconstructed whole-sky 
smoothed density field for a slice of 8000 km/sec wide, 
centered on the supergalactic plane, in z-space and in 
3D space, respectively. Note that the same Gaussian 
smoothing and cell units are used as in figures lb and 
Id. A different representation of the cluster density in 
this slice is also presented in Plate 1, where we plot a 
3D surface-plot visualization of the density field (similar 
to that presented by Dekel 1994), again in both z-space 
(upper panel) and 3D (lower panel). As also outlined 
by Scaramella (1994) who derived and discussed simi- 
lar maps but only in redshift space, the sparseness of 
the cluster distribution and the heavy smoothing could 
miss some features of the cosmic density field. Never- 



theless, our reconstruction technique should reveal the 
main large-scale features of the 3D cosmic density field. 

From these plots it is evident that eliminating the 
distortions due to peculiar velocities suppresses signif- 
icantly the amplitude of the density peaks in 3D with 
respect to that in the z-space. This effect is particularly 
evident in the Virgo-Hydra-Centaurus (or else Great At- 
tractor) region [centered at (X,Y) = (8.5,11)], in the Leo 
region [centered at (X,Y) = (11,14)] and around the 
Ursa-Major supercluster [centered at (X,Y)=(15,18)]. 
Interestingly, in the Perseus-Pegasus region [centered at 
(X,Y)=(13,10)] the opposite is true, ie., the 3D density 
amplitude is slightly higher than in the z-space case. 
This could be easily understood if there exists a coher- 
ent gradient of negative peculiar velocities, increasing 
in amplitude towards the near side of this structure (to- 
wards smaller X's, ie., towards the LG). 

Also, the shape of the density peaks, in the 3D case, 
appears to be slightly elongated along the line of sight. 
This effect, which is mostly prominent in the Great At- 
tractor and in the Shapley concentration [centered at 
(X,Y)=(3.5,14)], arises from the fact that in the linear 
regime the infall peculiar velocities within high density 
regions are coherent and tend to twist the isodensity 
contours along the line of sight. 

In figure 4 we present the 3D cluster density field 
for the four slices, each of 8000 km/sec width, above and 
below the supergalactic plane. Panel (a) and (6) corre- 
spond to the 4000 < Z < 12000 km/sec and 12000 < 
Z < 20000 km/sec slices respectively, while panel (c) 
and (d) represent the corresponding slices but for nega- 
tive Z's. Our intention is not to present a detailed anal- 
ysis of the structures evident in our contour plots (for 
such a task see TuUy 1987 and TuUy et al. 1992). How- 
ever, we would like to emphasize the fact that, most 
structures found in the supergalactic plane extend not 
only to the first 8000 km/sec wide slice, above or be- 
low the supergalactic plane, but also to the second such 
slice. For example the overdensity corresponding to the 
Shapley concentration extends to both negative Z slices. 
Also the Great Attractor region is connected to the Le- 
pus region [centered at (X,Y)=(8.5,10) in panel c] which 
seems to extend even further to the next slice (panel 
d) while a similar behaviour is found also for Perseus- 
Pegasus region. From panels (a) and (6) we see that 
the same is true for the Ursa Major supercluster, the 
Grus- Indus region [centered at (X,Y)=(3,2.5) in the su- 
pergalactic plane slice], the Leo region which seems to 
be connected with Hercules [centered at (X,Y)=(10,14) 
in the (a) and (6) slices] as well as the Pegasus-Pisces 
region [centered at (X,Y)=(17.5,4)] which is evident in 
the supergalactic plane and the two positive Z slices. 
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3.5 Redshift space distortions of the dipole 
vector 

One of the major scopes of this work is to compare the 
cluster dipole as measured in z-space with the recon- 
structed one in 3D-space. A simple eye inspection of 
figures 3 a,b and Plate 1 revealed that the peaks of the 
true density field have larger amplitudes when observed 
in redshift space. However, to see how these distortions 
affect the estimate of the cluster dipole require a quan- 
titative discussion. From a theoretical point of view the 
problem has already been addressed by Kaiser (1987). 
He showed that gravitationally driven peculiar veloci- 
ties can have a large effect on the peculiar acceleration 
vector of an object. If Qr is the amplitude of the true 
peculiar acceleration and Qz is the z-space peculiar ac- 
celeration, then the relationship between Qr and Qz de- 
pends both on the selection function and on the on the 
power spectrum considered. If one considers each power 
mode separately then, in the limit kRmax ^ 1 {Rmax 
is the depth of the survey and k is the wavenumber). 
Kaiser has shown that 



Hpc) 



(20) 



Qz - Qr = ffr— |— [1 + E{k)\ , 
where 

E{k) ~ \og[{kRma.f<t>{Rma.)l<t>{k~^)\ ■ 



(16) 



(17) 



Clearly the sign of E(k) depends on the selection func- 
tion and on the sample's depth. Although a precise E(k) 
evaluation require a good determination of the selec- 
tion function (that, for the clusters case, also depends 
on the relative weight), we can be confident that in our 
case E(k) > 0. In fact for the Abell and AGO cluster 
sample <t>(T) ~ 1 within the sample (the present analy- 
sis is limited to a depth of Rmax = 200 Mpc). In 
the linear regime the final value of Qz — Qr can be ob- 
tained by combining each mode separately (for large k 
the nonlinear evolution couple the different modes, and 
this approximation does not hold). Obviously to prop- 
erly combine the different modes requires the knowledge 
of the density fluctuation power spectrum. However, un- 
less one considers cosmological models with very large 
power on small scales, which are in conflict with current 
observations, we expect that for our cluster sample: 

ffz-ff. >0 (18) 



4 THE CLUSTER DIPOLE 

The dipole vector, for an observer placed at the center 
of the coordinate system, is defined as 



D = A/ 



Pc(r) — dr. 



(19) 



For our discrete cluster field density field, smoothed 
with a top hat window, the dipole becomes: 

3Mc \ - / Mi \ 



D 



47r 



■E 



(21) 



The monopole of the smoothed field is defined as: 

i 

Since the value of D/Af has a smaller scatter than that 
of D (see PV91) it can be used to compute the LG pe- 
culiar velocity, assuming that the dipole has converged 
to its asymptotic value within the volume sampled. In 
this case the relationship between the peculiar velocity 
and D/M is: 

, D(7-) HqR 

conv 



y{r)=p- 



(23) 



which, in linear theory is related to the peculiar velocity 
by 



3 M{Rconv) ' 

where Rconv is the final dipole's convergence depth. 
However, based on Montecarlo experiments, we found 
that a better estimate of the peculiar velocity of the 
observer is obtained directly from eq.(13). The average 
1(7 scatter between v, computed from eq.(23), and the 
true one is ~ 12% if 850 'clusters' are used, a number 
comparable to that of our cluster (true and synthetic) 
sample. The discrepancy originates from the fact that 
the monopole is a good estimator of the average density 
only in the limit of continuous density field. Using 8500 
objects, for example, reduces the scatter to 7%. This 
effect could introduce a further uncertainty in the LG 
peculiar velocity computed via eq.(23). 

Gomputing the dipole and the monopole separately 
is nevertheless important for assessing the effective con- 
vergence of the peculiar velocity vector since a neces- 
sary condition for final convergence is that while the 
dipole converges, the cumulative monopole continues to 
grow linearly with distance. This condition, however, is 
not sufficient to guarantee the final convergence of the 
dipole, which could still increase after having reached a 
plateau (cf. Plionis, Goles &; Gatelan 1993). 

4.1 Possible Systematic Effects 

The reconstructed cluster distribution could depend on 
a number of free parameters that are only partially de- 
termined by observational constraints and theoretical 
requirements. We systematically explored the influence 
of the parameter choice by allowing each free param- 
eter to vary within an experimentally plausible range. 
In practice the reconstructed cluster dipole, and there- 
fore the predicted LG peculiar velocity, can be regarded 
as a function of many variables (the free parameters) 
that are listed in table 1. Boldface quantities refer to 
what we call the standard case, defined by the following 
choice of parameters: (a, a) = (0.3,0.2), = 20 
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Mpc, Qo = 0.5, (3 = 0.25, direct homogenization scheme 
(eq.5) with the 200 Mpc sphere being divided into 
10 equal volume shells. For each of the 48 possible com- 
binations of the free parameters, to which we refer to as 
"models", we performed 10 independent reconstructions 
and we computed the corresponding 3D dipole. To ex- 
plore the influence of each variable on the predicted LG 
dipole separately we adopted the following procedure: 

• The parameter of interest is selected and models 
are divided according to the value chosen for that 
parameter (e.g. qo = 0.2 or 0.5). Models having 
the same parameter values are grouped in the same 
category. 

• The average dipole at 170 Mpc is computed for 
all the models in the same category and the same 
procedure is repeated for all different categories. 

• The different average dipoles are compared. The 
differences in amplitude and the misalignment an- 
gles between different average dipoles are listed in 
table 2. 

It is evident from table 2 that the only appreciable 
systematic effect comes from the choice of the (3 param- 
eter, as we will discuss below. Varying the other exper- 
imental parameters (i.e. the galactic absorption coeffi- 
cients, the smoothing radius and the homogenization 
procedure) generates systematic errors much smaller 
than the intrinsic uncertainties of the reconstruction. 

4.1.1 The effect of varying the (3 parameter 

To appreciate the influence of the (3 parameter on the 
reconstruction process we have to disentangle two dif- 
ferent effects: the role of the deceleration parameter, go, 
that enters through eq.(l) in defining the outer limit of 
the sample, and that of the (3 factor, whose arbitrary 
initial value has to be specified for the LIRM to be ap- 
plied. The first effect turned out to be much smaller 
than the second one. 

The effect of varying the (3 parameter in the LIRM 
was explored by performing the reconstruction with the 
same qo value but using two different (3 values (0.25 
and 0.114). Although these values are different by a 
factor ~ 2, the asymptotic amplitudes of the recon- 
structed dipole are very similar, differing only by 8%, 
being smaller for the standard (3 = 0.25 case, while 
the direction does not change appreciably. The complete 
analysis was performed by exploring two different class 
of models having (/3 = 0.25, qo = 0.5) and (/3 = 0.114, 
qo = 0.2), respectively, confirmed the above results (see 
table 2). It is therefore clear that there is a small, al- 
though systematic, effect related to the value of the (3 
parameter which enters as an input in the reconstruc- 
tion algorithm. 

It is interesting to compare this result to that ob- 
tained by Strauss et al. (1992, [ hereafter SYDIIF92]). 



The asymptotic amplitude of their IRAS galaxy dipole, 
on scales larger than ~ 100 h~^ Mpc, strongly depends 
on the input (3 value. The reason for this behaviour was 
ascribed to redshift space distortions, to which they re- 
fer as the 'Kaiser effect', that in their case is mainly 
due to the LG motion with respect to the distribution 
of distant galaxies. In this simple situation the prob- 
lem is similar to that of a 'rocket-born' observer who 
resides in a uniform Universe at rest with respect to 
the CMB on large-scales and whose peculiar accelera- 
tion originates locally (Kaiser &; Lahav 1988). In this 
idealized case the redshift space distortions apprecia- 
bly affect the measured peculiar acceleration only if the 
selection function falls below one, well within the sam- 
ple's depth. This is the case for the IRAS galaxy sample 
studied by SYDIIF92, where < 1 well within 100 h~^ 
Mpc. For our sample, however, P{z) is close to unity 
within 200 h~^ Mpc so that the 'Kaiser Effect' is much 
less important, as demonstrated by the stability of our 
dipole against the variations of (3; especially at large 
depths where the sensitivity to the input /3-parameter 
would be manifested by variations of the dipole ampli- 
tude shape as a function of (3 (cf. Fig. 6 of SYDIIF92). 

4.2 Error Estimate 

To obtain a reliable error estimate of the reconstructed 
dipole we need to account for several different sources of 
errors. We estimate the intrinsic error, ai, in the recon- 
struction algorithm from the 10 independent reconstruc- 
tions performed for each model explored. Other sources 
of errors, which will be described below, are the obser- 
vational error (Jo, the 'shot-noise' error (Js„ and the zero 
point uncertainty (Jzp. In what follows we will assume 
that all these errors are independent so that the total 
uncertainty ctt will be simply estimated by adding the 
errors in quadrature. 

4.2.1 Observational error 

We define as observational error the uncertainty derived 
from the different free parameter choice, listed in table 
1, whose separate influence on the reconstructed dipole 
was explored in the previous subsection. The recon- 
structed dipole can be therefore regarded as a function 
of the smoothing length, the homogenization scheme 
and its binning, parameters that we allowed to vary in 
a plausible range of values. The galactic absorption set 
and the (3 factor, instead, are considered fixed by obser- 
vations to (0.2,0.2) and to (3 = 0.25, respectively. We 
chose (3 = 0.25 since, as we will see in section 5, this 
value is closer to that one derived from the comparison 
with the CMB dipole. We found that varying these pa- 
rameters does not produce any systematic effect in the 
dipole determination and thus we can consider them 
as Gaussian variables and the dipole as a multivariate 
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Gaussian function of the explored parameters. Follow- 
ing this argument we compute (Jo as the variance of 
the dipole obtained by averaging over all models having 
(a, a) = (0.3,0.2) and (3 = 0.25 (our standard case). 
The observational error at 170 Mpc produces an 

uncertainty of ~ ±65 km/sec in the dipole am- 

plitude and of ±4° and ±2° along b and I directions 
respectively. 

4-2.2 Shot-Noise error 

The usual shot-noise estimate assumes that luminous 
objects have been drawn by a Poisson process from 
the underlying density field. However the reconstruc- 
tion procedure relies on the hypothesis that clusters of 
galaxies are biased tracers the mass. A more meaning- 
ful and self-consistent estimate of the shot noise error 
has been proposed by SYDHF92 that accounts for two 
different effects: 

• The Poissonian error due to the fact that luminous 
objects become random tracers of the density field 
when the selection function falls below unity. 

• An uncertainty due to the fact that cluster masses 
are randomly taken from an underlying mass dis- 
tribution. In our case we estimated the mass of real 
clusters from their Abell listed number of galaxies 
while for synthetic clusters, which constitute ~ 50% 
of the whole cluster population, the mass was set 
equal to the average real cluster mass. 

Using our formalism, the SYDFF92 shot noise com- 
puted for the Wrei = 1 case, reduces to: 

i 

where the sum extends only over the Ns synthetic clus- 
ters. Tic is the average cluster density defined from 
eq.(14) and K = — 1, with fii being the mass of 

the real clusters in units of average mass. From the real 
clusters mass variance we have K=0.25. To compute the 
influence of the shot noise on the amplitude and the 
direction of the dipole we will assume that each com- 
ponent of the dipole is a Gaussian with zero mean and 
o'sn,3: = o'sn,y = o'sn,z- The l-dimcnsioual shot noise er- 
ror at 170 Mpc reflects m an ~ ±60 /3"^ km/ sec 
amplitude and a ~ ±4° directional uncertainty. 

4-2.3 Zero point error 

The gravitational influence on the LG motion from 
nearby clusters, such as the Virgo cluster which was not 
included in the Abell/ AGO catalogue due to its proxim- 
ity and therefore of its large angular dimension, are not 
negligible. Neglecting them in the cluster-dipole analy- 
sis, causes a zero point offset on the predicted LG pecu- 
liar acceleration. A possible way to overcome this prob- 



lem is to merge the cluster dipole with that of some 
galaxy sample that trace more densely the local mass 
distribution (Scaramella Vettolani &; Zamorani 1994, 
[SVZ94 hereafter]). With this approach, however, the 
offset can be evaluated only assuming a priori a value 
for the relative bias of the cluster and galaxy popula- 
tions, whose value is very uncertain (for a complemen- 
tary approach see Plionis 1995 in preparation) - 

Here we adopt a simpler phenomenological ap- 
proach by identifying the zero point offset with a 200 
km/sec LG Virgocentric infall to which we attach a 
±100 km/sec uncertainty. This results in an uncertainty 
of ^'j''^ degrees and ^^j^ degrees in the dipole direction 
along I and b, respectively. Analogously, for the dipole's 
amplitude we obtain an uncertainty of km/sec. 

4.3 The 3D Cluster Dipole 

After reconstructing the cluster distribution with the 
RSCM and the LIRM, we computed the cluster dipole 
both in redshift and real-space. Figure 5a shows the am- 
plitude of the peculiar velocity for the standard model, 
computed using eq.(13), as a function of distance. Open 
dots refer to the z-space dipole while filled symbols 
represent the reconstructed 3D cluster dipole. We also 
plotted, for reference, the dipole computed without any 
Abell/ACO homogenizing scheme (i.e. Wrei = 1) in z- 
space (continuous line). Errorbars represent 1 a total 
errors. 

Our main result is that the asymptotic value of the 
3D dipole is significantly less than the z-space dipole, 
as predicted by eq.(18). Removing redshift space dis- 
tortions erases the artificial redshift space clustering, 
leading to a smaller, by ~ 23%, asymptotic dipole am- 
plitude. Although the amplitude of the 3D dipole is sig- 
nificantly less than what previously found in cluster- 
dipole z-space studies (PV91; SVZ91), the qualitative 
dipole behaviour is similar in z-space and 3D-space: a 
"bump" of the velocity amplitude around 50h~^ Mpc 
followed by a decrease due to the competing pull of 
the Great Attractor with the Perseus-Pegasus regions, 
a secondary increase after a plateau and an asymptotic 
convergence beyond ~ 170 Mpc once that the Shap- 
ley concentration has entered into the sampled volume. 
The mass distribution beyond 60 Mpc is respon- 

sible for ~ 30% of the total LG peculiar acceleration. 
As already noticed from the qualitative analysis of the 
density field maps, the effect of eliminating the appar- 
ent z-space clustering is particularly important in lo- 
cal high density regions, such as the GA area, which 
is clearly responsible for the significant decrease of the 
~ 50 Mpc bump in the 3D cluster dipole case. In 
Figure 5a we also plot, in arbitrary units, the 3D space 
monopole term, which grows linearly out to the depths 
sampled which indicates that the dipole convergence at 
~ 170 Mpc is not an artifact of ill sampling. 
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Figure 5b displays the cumulative direction of the 
LG velocity. The starred symbol indicates the CMB 
dipole apex, while the skeletal symbol represents the 
CMB direction after correcting the LG for a Virgocen- 
tric infall of 200 km/sec. Errorbars represent the 1 a 
total errors. The relatively large error (~ 10°) in the 
dipole direction, along b, is due to the uncertainty re- 
lated to the amplitude of the Virgocentric infall. The 
reconstructed 3D dipole points ~ 10° away from the 
CMB apex when the Virgo infall is taken into account. 
No significant differences between the 3D and z-space 
dipole directions was found. 



4.3.1 The effect of varying the Mask model 

It is interesting to investigate the influence of the Mask 
model, adopted to reconstruct the whole sky z-space 
cluster distribution, on the reconstructed dipole. Here 
we compare the RSCM reconstructed standard dipole 
with the two obtained by using the same set of stan- 
dard parameters but adopting the RHCM and RTCM 
mask models. In both cases there is a small increase of 
the dipole's amplitude, by ~ 3% — 4%, and the direction 
points ~ 4 ± 4° away from the original one (uncertain- 
ties refer to the intrinsic errors). We conclude that no 
significant dependance was found on the Mask choice, 
which is to be expected since all of them are based on 
a cloning scheme. 

A more relevant comparison is probably between 
two conceptually different methods. For example be- 
tween our RSCM and the spherical harmonic recon- 
struction method (cf. Lahav 1987; Yahil et al 1986; Plio- 
nis 1988 and PV91 for use of this method in the cluster 
dipole context). In Figure 6 we compare our standard 
RSCM ^-reconstructed dipole with that derived using 
the PV91 method for the same sample and the same set 
of parameters. In the insert we plot the relative velocity 
fluctuations, defined as 2 x {vrscm — vsh)I{vrscm + 
vsh). It is evident that both methods give identical re- 
sults, except in the inner volume 50 h~^ Mpc) where 
the number of clusters is extremely small and therefore 
shot noise effects are large. This is quite extraordinary 
given the completely different methods used. 

Finally, we explored the effect of using the MSS 
with which the reconstruction scheme adopted in the 
inner region 200 h~^ Mpc) is extended up to the 
limiting sample depth (250 h~^ Mpc). With the MSS 
the 3D dipole differs from the standard one only be- 
yond 200 h~^ Mpc, where the amplitude monotonically 
increases (by ~ 8%) and the direction systematically 
drifts towards lower galactic latitudes. As we discussed 
in section (3.1), using the MS'S' increases the shot noise 
errors in a region where the selection function is signif- 
icantly smaller than 1. The shot noise errors are proba- 
bly responsible for the drift in the dipole direction, while 



they also enhance the 'Kaiser Effect' which causes the 
observed amplitude increase. 



5 COSMOLOGY FROM THE CLUSTER 
DIPOLE 

5.1 Constraints on the value of (3 

From the measured asymptotic value of the cluster 
dipole it is possible to constrain the value of /3. An 
estimate of this parameter can be obtained by simply 
assuming that the cluster 3D dipole has converged at 
~ 170 h~^ Mpc, where it flattens, in accordance with 
SVZ91 and PV91. Using eq.(13) to compare the dipole 
amplitude with the LG motion, corrected for the Virgo- 
centric infall leads to a value of the /3-parameter: 

/3 = 0.2liroI 

If we erroneously did not account for the Virgo infall 
we would obtain /3 = 0.25^q qs- Errors come from total 
dipole uncertainties, as discussed in section 4. 

Although the dipole convergence at ~ 170 h-^ Mpc 
is real, since at that distance the cluster monopole is 
still linearly increasing, it is impossible to prove that the 
dipole has converged to its final value. The only way out 
is to evaluate the amplitude of the external contribution 
to the observed dipole using statistical methods. 

Juskiewicz, Vittorio &; Wyse (1990) computed the 
probability distribution of the peculiar velocity of a 
generic observer, constrained to be moving with a ve- 
locity of 620 km/sec with respect to the CMB, as a 
function of the sampled volume. From this distribution 
they computed the expectation value for the estimated 
LG velocity and its variance (their eq. 8a and 8b). The 
expectation value at 170 h~^ Mpc weakly depends on 
the cosmological model adopted. We explored a number 
of possible models ranging from the Standard CDM one 
to the phenomenological model proposed by Branchini 
et al. (1994) that resembles a tilted CDM with n = 0.7 
and to the Peacock &; Nicholson (1994) spectrum with 
more power on large scale (similar to a standard CHDM 
model). In all these cases the expectation value for the 
LG velocity, at 170 h~^ Mpc coincides with the true one. 
Also the variance slightly depends on the cosmological 
model, being larger in correspondence to the Peacock &; 
Nicholson (1994) spectrum (~ 600 /3"^ km/ sec that, to 
be conservative, will be used as the typical cosmological 
variance). A probably better, although exaggerated, es- 
timate of the (3 uncertainty can be therefore obtained by 
adding in quadrature the cosmological variance to the 
total error. The resulting (3 value is thus constrained 
into the range 

The 1 a range is wide enough to overlap with previous 
estimate of (3 based on the z-space cluster distribution 



Reconstructing Positions & Velocities of Clusters I. 13 



(PV91, SVZ91, SVZ94). Assuming fio = 1 we obtain a 
value for the cluster bias parameter 

According to the linear biasing prescriptions if Ba 
and 6b are the biasing factors of objects A and B with 
respect to the matter with 6a > 6b, the following re- 
lation holds: Bab ~ bA/bB, where 6ab is the biasing 
factor of objects A with respect to objects B. If the bi- 
asing factor of Abell/ACO clusters with respect to IRAS 
galaxies is bd ~ 4.5 (cf. Peacock &; Dodds 1994; Jing 
&; Valdarnini 1993) we have that the above 6c value is 
consistent with a bias parameter of IRAS galaxies with 
respect to the mass of 6/ ~ I.I^q g- 

5.2 Constraining Dark Matter Models 

The observed cumulative acceleration acting on the LG 
can be used to discriminate between different cosmolog- 
ical models. Useful constraints can be imposed by re- 
quiring a cosmological model to reproduce the observed 
cumulative dipole. As outlined by SVZ94, the dipole in- 
crease in the range beyond 80 Mpc is a characteris- 
tic feature which can be used as a cosmological test. To 
quantify the constraint we consider the ratio of peculiar 
velocity X(Ri,R2) induced on an observer placed on the 
origin from fluctuations within a sphere of radius Ri to 
that from within a sphere of radius R2 > iii . Juskiewicz, 
Vittorio &; Wyse (1990) computed the probability dis- 
tribution for this ratio as a function of the cosmologi- 
cal model through the power spectrum P{k). Here we 
computed the probability density distribution for ratio 
X(Ri,R2) °f ^he peculiar velocities generated within the 
two spheres of Ri = 100 and R2 = 160 Mpc. The 
probability distributions computed for the three cosmo- 
logical models previously used, are displayed in Figure 
7. The continuous line represents the experimental val- 
ues from the 3D standard dipole in Figure 6 assuming 
a Gaussian distributed total errors. Short dashed line 
represents the CDM standard model, long dashed line 
refers to the Branchini, Guzzo &; Valdarnini (1994) phe- 
nomenological model while the Peacock and Nicholson 
(1991) model is displayed as a dot-dashed line. As al- 
ready noticed by SVZ94, models with a large scale of 
coherence are preferred by the observed cluster dipole 
even when redshift space distortions are removed. A de- 
tailed analysis, using large Zel'dovich simulations of the 
cluster distribution, on the constraints put by the clus- 
ter dipole on the different Dark Matter models will be 
presented in Tini Brunozzi et al. (1995). 

6 CONCLUSIONS 

In this work we described a self consistent method to re- 
construct the 3D distribution and peculiar velocities of 



galaxy clusters from their observed positions in the red- 
shift space. Our approach is valid within the framework 
of gravitational instability and relies on the assump- 
tion of linear theory and linear biasing. The technique 
presented aims at recovering, in a statistical sense, the 
distribution of unobserved clusters and corrects the ap- 
parent cluster positions for peculiar velocity driven red- 
shift space distortions. Tests performed using synthetic 
cluster distributions, obtained from numerical experi- 
ments, showed that the intrinsic reliability of this tech- 
nique allows the true cluster positions to be recovered 
with an average error of ~ 1.6 %, much better than 
the best available experimental determination of cluster 
distances (> 10%). H owever, the actual reconstruction 
uncertainty is larger since we adopted a phenomeno- 
logical method, to reconstruct the whole sky z-space 
distribution of clusters, characterized by several param- 
eters with associated observational errors. As necessary 
in any phenomenological approach, we tested the sta- 
bility of the final results on the various experimental 
parameters and found them to be surprisingly robust. 

Our main result is that redshift space distortions 
cause an artificial increase of the cluster dipole ampli- 
tude by ~ 23%, while they do not affect appreciably the 
dipole direction which points ~ 10° away from the CMB 
apex. The cumulative 3D dipole qualitative matches the 
behaviour of the z-space one: there is a bump at ~ 50 

Mpc and, after a plateau, it rises by ~ 30% and 
stabilizes beyond 170 Mpc where it reaches a real, 
although not necessary final, convergence. The reduc- 
tion of the 3D dipole amplitude with respect to that of 
the z-space one arises mainly from the removal of the 
strong apparent clustering in the Great Attractor and 
Shapley regions and its enhancement in the Perseus- 
Pegasus region, as can be clearly seen in the smoothed 
density maps of Figure 3a,b and Plate 1. These results 
are robust in the sense that a search for the presence 
of systematic biases gave negative results. In practice 
we found that the freedom of modelling the galactic ab- 
sorption, the radial selection function, the ZoA filling 
scheme, the homogenization and the smoothing proce- 
dures cause systematic errors much smaller than the 
intrinsic error in the reconstruction procedure. Interest- 
ingly, since the selection function of our cluster sample 
is ~ 1 up to ~ 190 Mpc, our 3D dipole determi- 
nation is only marginally affected by the 'Kaiser effect' 
which has been found to be a serious problem in the 3D 
galaxy dipole estimations. 

From the asymptotic dipole amplitude it is possible, 
in principle, to determine the (3 parameter. The intrinsic 
cosmological variance, however, reflects in a large uncer- 
tainty whose magnitude depends on the unknown cos- 
mological background. In the conservative and widely 
accepted hypothesis that the cosmic density field is char- 
acterized by a large correlation length, we estimated 
/3 = 0.21^Q 'l'^ . This value is larger than previously found 
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based on the analysis of the z-space cluster distribution, 
the reason being that we have accounted for redshift 
space distortions. The large la uncertainty, however, 
does not allow us to set stringent limits on the value of 
(3. Better constraints will be presented in a forthcom- 
ing paper in which we will analyze the reconstructed 
cluster peculiar velocities and compare them to those 
determined observationally. 
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Table Captions 

Table 1: Experimental parameters. For each parameter 
the values adopted are shown. Each model explored is 
characterized by a different set of parameters. Boldface 
quantities represent the 'Standard Model' set. 

Table 2: Differences between the different average 
dipoles. First column displays the parameters defining 
the two compared categories. Column 2 shows the per 
cent difference of the dipole amplitude at 170 Mpc 
together with the relative standard deviation. Columns 
3 and 4 contain the average misalignment and its stan- 
dard deviation at 170 Mpc along I and b, respec- 
tively. 



Figure Captions 
Figure 1: 

(a) and (c): The projected whole-sky z-space cluster 
distribution in Cartesian supergalactic coordinates. 
Filled symbols represent real Abell/ACO clusters 
while open symbols the synthetic objects. Short 
dashed lines delineate the |6| < 20° region, (a) 
Xs„p-Ys„p projection and (c) Xs„p-Zs„p projection. 

(b) and (d): The projected smoothed whole-sky z- 
space cluster density field (of panels (a) and (c)). 
A 2-D Gaussian kernel was used with Rsm = 2000 
km/sec. (b) Xs„p-Ys„p projection and (d) Xs„p- 
Zsup projection. 

Note that cell units are used (see text for definitions). 



Figure 2: The distribution of relative percentage dif- 
ferences, £i, between the true and reconstructed posi- 
tions of simulated clusters (see text for details) which 
were used to assess the reliability of our reconstruction 
method. 



Figure 3: The projection onto the supergalactic plane 
of the smoothed cluster density field with —4000 < 
Zsup < 4000 km/sec. (a) z-space (b) reconstructed 3D- 
space. Same smoothing procedure used as in figure 2. 

Figure 4: Xs„p-Ys„p projections of the smoothed 3D- 
space cluster density field for 4 slices in Zs„p. (a) 4000 < 
Zs„p < 12000 km/sec. (b) 12000 < Zs„p < 20000 
km/sec. (c) -12000 < Zs„p < -4000 km/sec. (d) 
-20000 < Zs„p < -12000 km/sec. 

Figure 5: Cluster dipole. 

(a) Dipole amplitude as a function of distance with 
open symbols representing the reconstructed whole- 
sky z-space dipole and filled ones the reconstructed 
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3D dipole. The continuous thin line refers to the z- 
space dipole without any Abell/ACO homogenizing 
scheme (i.e. Wrei = 1). In arbitrary units, we plot 
also the 3D-space monopole (small filled squares), 
(b) Dipole misalignment angle. The starred symbol in- 
dicates the CMB dipole apex, while the skeletal 
symbol represents the same but after correcting for 
a Virgocentric infall component of the LG motion, 
with an amplitude of 200 km/sec. 

Errorbars are 1 a total uncertainty (see text for details). 

Figure 6: Comparison of the z-space dipole derived by 
using our standard RSCM method to reconstruct the 
whole-sky cluster distribution and that derived using 
the spherical harmonic method (cf. PV91) for the same 
sample and the same set of parameters. In the insert we 
plot the relative velocity fluctuations between the two 
determinations. 

Figure 7: The probability density distribution for the 
dipole amplitude ratio between the two spheres with 
radii Ri = 100 and R2 = 160 Mpc. The contin- 

uous line represents the experimental values from the 
3D standard dipole of Figure 7. Short dashed line rep- 
resents the CDM model, long dashed line refers to the 
Branchini, Guzzo &; Valdarnini (1994) phenomenologi- 
cal model while the Peacock and Nicholson (1991) model 
is displayed as a dot-dashed line. 



PLATE 1: Surface-plot visualization of the density 
field of the slice -4000 < Zsup < 4000 km/sec; the up- 
per panel represents the density field in z-space and the 
lower panel the corresponding 3D-space density field. 
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